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Abstract 

Objectives: This meta-analysis was conducted to investigate whether the status of rheumatoid factor (RF) and anti-cyclic 
citrullinated peptide (anti-CCP) antibody are associated with the clinical response to anti-tumor necrosis factor (TNF) alpha 
treatment in rheumatoid arthritis (RA). 

Methods: A systemic literature review was performed using the MEDLINE, SCOPUS, Cochrane Library, ISI Web of Knowledge, 
and Clinical Trials Register databases, and Hayden's criteria of quality assessment for prognostic studies were used to 
evaluate all of the studies. The correlation between the RF and anti-CCP antibody status with the treatment effect of anti- 
TNFot agents was analyzed separately using the Mantel Haenszel method. A fixed-effects model was used when there was 
no significant heterogeneity; otherwise, a random-effects model was applied. Publication bias was assessed using Egger's 
linear regression and a funnel plot. 

Results: A total of 14 studies involving 5561 RA patients meeting the inclusion criteria were included. The overall analysis 
showed that the pooled relative risk for the predictive effects of the RF and anti-CCP antibody status on patient response to 
anti-TNFa agents was 0.98 (95% CI: 0.91-1.05, p = 0.54) and 0.88 (95% CI: 0.76-1.03, p = 0.1 1), respectively, with I 2 values of 
43% (p = 0.05) and 67% (p<0.01), respectively. Subgroup analyses of different anti-TNFa treatments (infliximab vs. 
etanercept vs. adalimumab vs. golimumab), response criteria (DAS28 vs. ACR20 vs. EULAR response), follow-up period (>6 
vs. <6 months), and ethnic group did not reveal a significant association for the status of RF and anti-CCP. 

Conclusions: Neither the RF nor anti-CCP antibody status in RA patients is associated with a clinical response to anti-TNFa 
treatment. 
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Introduction 

Rheumatoid arthritis (RA) is a chronic inflammatory autoim- 
mune disease that affects approximately 1% of the population 
worldwide [1]. OLE_LINK10Although the introduction of anti- 
TNFa agents has dramatically improved the outcome of RA, there 
unfortunately remains a proportion of RA patients who do not 
exhibit an adequate response to this treatment. Considering the 
high cost and potential side effects of anti-TNFa treatment, it is 
important to identify those RA patients who will be more likely to 
respond to these agents. Indeed, numerous studies have been 



conducted to investigate potential predictors for patient response 
to anti-TNFa therapy [2-4] . 

Both rheumatoid factor (RF) and antibodies against cyclic 
citrullinated peptide (anti-CCP) are regarded as serological 
markers of RA [5,6] . Some studies have suggested that the status 
of RF or anti-CCP antibody in RA patients is associated with a 
clinical response to anti-TNFa treatment [7-14], whereas such a 
correlation was not found in other studies [15-19]. Thus, no 
definite conclusion has been reached to date. 
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Flowchart of the selection process of included studies for meta-analysis 
Bolded risk ratio highlight statisti callysignificant results. 



1626 records were identified 
through database searching 



23 additional records were identified 
using "search alert" function 



1649 records were identified through 
systematic literature search 



1296 records screened on 
basis of title and abstract 



353 records excluded 
as duplicates 



1225 records excluded with 
reasons: 

• Patients, interventions, 
outcomes didn't meet 
general inclusion criteria. 

• Serotype of RF or anti-CCP 
of patients wasn't reported. 



71 records assessed on 
basis of full-text 







15 articles were included 
in qualitative synthesis 







14 Studies included in 
qualitative synthesis 
(meta-analysis) (total 
numberof patients = 
5561) 



56 records excluded with 
reasons: 

• Patients, interventions, 
outcomes, and data 
reported didn't meet 
inclusion criteria(48) 

• Conference paper(l) 

• Review article(6) 

• Lettertotheeditor(l) 



Figure 1. Flowchart of the study selection. 

doi:1 0.1 371 /journal.pone.0089442.g001 

We performed a meta-analysis to investigate whether RF and 
anti-CCP have predictive value for a clinical response to anti-TNFa 
treatment. Suitable studies investigating an association of the status 
of RF or anti-CCP and response to anti-TNFa treatment were 
searched and included. We also performed subgroup analyses on 
different variables to explore potential sources of independent 
predictive factors for an effect of anti-TNFa treatment. 

Methods 

Search strategy 

A literature search was performed for all studies evaluating an 
association between the status of RF or anti-CCP antibody and a 



response to anti-TNFa therapy in RA patients using the Medline, 
Cochrane Library, SCOPUS (including EMbase), ISI Web of 
Knowledge, and Clinical Trials Register (clinical trials.gov) 
databases. The following keywords were searched: rheumatoid 
arthritis, anti-TNFa, rheumatoid factor, anti-cyclic citrullinated 
peptide antibody, clinical trials, and systematic review. Synonyms 
and spelling variations were taken into account (Search strategy for 
Scopus was listed in Table SI in File SI). There was a limitation 
with regard to language, i.e., we only considered English 
publications, but not the year of publication. We also contacted 
authors to request a full-text review or specific data from studies 
when there was no electronic version of the full text or sufficient 
data for the meta-analysis. Citations were reviewed to search 
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Study or Subgroup 



Risk Ratio 
M-H. Random. 95% CI 



Risk Ratio 
M-H. Random. 95% CI 



Alexandras2009 

Bobbio2006 

Canha~o2012 

Hyrich2006 

Keystone2009 

Lequerre2006 

Lilian2011 

Mancarella2007 

Potter2009 

Vasilopoulos201 1 

Wijbrandts2008 

Wouter2008 



0.99 [0.58, 1.70] 
0.99 [0.71, 1.37] 
1.03 [0.86, 1.23] 
1.02 [0.97, 1.07] 
1.54 [0.93,2.55] 
0.92 [0.65, 1.31] 
0.89 [0.76, 1.06] 
0.73 [0.58, 0.90] 
0.94 [0.83, 1.07] 
0.88 [0.66, 1.17] 
1.48 [0.99, 2.22] 
1.06 [0.91, 1.23] 

0.98 [0.91, 1.05] 



Total (95% CI) 
Total events 

Heterogeneity: Tau 2 = 0.01; Chi 2 = 19.36, df = 11 (P = 0.05); I 2 = 43% 
Test for overall effect: Z = 0.61 (P = 0.54) 



0.5 0.7 1 1.5 2 

Favours serotype negative Favours serotype positive 



R 












Risk Ratio 


Risk 


Ratio 


Study or Subqroup 




M-H. Random. 95% CI 


M-H. Random. 95% CI 


1.1.1 EULAR criteria 




0.99 [0.58, 1.70] 






Alexandras2009 








Bobbio2006 




0.99 [0.71, 1.37] 






Canha"o2012 




1.03 [0.86, 1.23] 






Hyrich2006 




1.02 [0.97, 1.07] 






Lequerre2006 




0.92 [0.65, 1.31] 






Mancarella2007 




0.73 [0.58, 0.90] 






Potter2009 




0.94 [0.83, 1.07] 


— » 




Wouter2008 




1.06 [0.91, 1.23] 






Subtotal (95% CI) 




0.98 [0.91, 1.05] 






Total events 










Heterogeneity: Tau 2 = 


0 


00; Chi 2 = 10.61, df = 7 (P = 0.16); I 2 = 34% 






Test for overall effect: 


z 


= 0.66 (P = 0.51) 






1.1.2 DAS28 criteria 










Lilian2011 




0.89 [0.76, 1.06] 


■ 




Vasilopoulos201 1 




0.88 [0.66, 1.17] 


■ — 




Wijbrandts2008 




1 .48 [0.99, 2.22] 




■ ■ 


Subtotal (95% CI) 
Total events 




1.02 [0.73, 1.41] 






Heterogeneity: Tau 2 = 


0 


06; Chi 2 = 8.63, df = 2 (P = 0.01 ); I 2 = 77% 






Test for overall effect: 


z 


= 0.10 (P = 0.92) 






1.1.3 ACR20 










Keystone2009 




1.54 [0.93, 2.55] 




■ 


Lilian2011 




0.95 [0.64, 1.39] 


■ 




Subtotal (95% CI) 
Total events 




1.18 [0.69, 2.01] 






Heterogeneity: Tau 2 = 


0 


10; Chi 2 = 2.84, df = 1 (P = 0.09); I 2 = 65% 






Test for overall effect: 


z 


= 0.61 (P = 0.54) 
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Figure 2. Forest graphs of the meta-analysis of RF status and response to anti-TNFo agents. The overall analysis of RF status showed a 
pooled RR of 0.98 (95% CI: 0.91-1.05, p = 0.54) and an I 2 of 43%. Subgroup analyses on different response criteria revealed no significant differences. 
doi:1 0.1 371 /journal.pone.0089442.g002 

relevant original studies, and an electronic search alert was set to Study Selection 

cover recent studies. There were 1649 references identified by the literature search. 

Three individual investigators (QL, YY, & XL) evaluated the 
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Table 3. Subgroup meta-analysis of RF and RA patient response according to different anti-TNFa agents, follow-up periods, 
response criteria, and ethnic groups. 













Test of association 




Test of heterogeneity 




Subgroup 


Population 


No. of 
studies 


No. of 
patients 


RR 


95%CI 


P value 


Model 


P value 


I 2 




overall 


12 


5374 


0.98 


0.91-1.05 


0.54 


R 


0.05 


43% 


Anti-TNFa agent* 


infliximab 


4 


1827 


1.03 


0.89-1.19 


0.71 


R 


0.28 


23% 




adalimumab 


2 


248 


0.98 


0.81-1.18 


0.8 


R 


0.09 


65% 




etanercept 


1 


1267 


1.05 


0.97-1.12 


0.21 


NA 


NA 


NA 




golimumab 


1 


178 


1.54 


0.93-2.55 


0.09 


NA 


NA 


NA 


Follow up period 


^6 months 


9 


5017 


0.96 


0.90-1.03 


0.27 


R 


0.13 


37% 




<6 months 


3 


357 


1.25 


0.87-1.78 


0.22 


R 


0.11 


54% 


Ethnic group 


European 


10 


5146 


0.98 


0.91-1.06 


0.59 


R 


0.10 


39% 




South American 


1 


50 


0.89 


0.76-1.06 


0.19 


NA 


NA 


NA 



R, random-effects model; NA, not applicable. 

*The number of studies and number of patients receiving individual anti-TNFa agent in this table are those that we could identified after literature review and 

contacting the authors. 

doi:1 0.1 371 /journal.pone.0089442.t003 



references, and the decision of inclusion was made by consensus. A 
study was included based on the following criteria: 1) the patients 
were older than 16 years old, diagnosed with RA using ACR 
criteria, and treated with at least one anti-TNFa agent 
(adalimumab, infliximab, etanercept, certolizumab, or golimu- 
mab); 2) efficacy was measured using EULAR or ACR or DAS28 
criteria after a minimum duration of 1 2 weeks; and 3) the status of 
RF or anti-CCP antibody at baseline and sufficient data to 
calculate the risk ratio (RR) were reported in the study. The 
following information was extracted from each study: the study 
design, patients' characteristics, baseline status of RF or anti-CCP 
antibody, interventions, outcomes, and study duration. 

Study assessment 

All of the studies included were evaluated for potential bias 
using Hayden's criteria of the quality assessment of prognostic 
studies [20] (The details of the modified criteria were provided in 
Table S2 in File S2). Six domains were considered in the bias 
assessment method: "study participation", "study attrition", 
"prognostic factor measurement", "outcome measurement", 
"confounding measurement and account", and "analysis". Scores 
of "0", "1", or "2" represented "high bias", "partial bias", and 
"low bias", respectively. If the available information was 
insufficient to obtain a decision after contacting the authors, 
"unsure" was used, and no score was calculated. The total scores 
were calculated to show the level of study quality; studies with a 
score of "1 1-12" were rated as having relatively high quality, with 
"9-10" as moderate, and "less than 8" as low quality. The 
evaluation was performed independently by three investigators 
(QL, YY, & XL) to improve the validity of the results. 

Statistical analysis 

We estimated the risk ratio of the RF or anti-CCP antibody 
status for a response to anti-TNFa treatment. The point estimate 
of the risk ratio (RR) and 95% confidence intervals (CI) were 
calculated for each study, and a Mantel Haenszel analysis was 
used to calculate the pooled RR and 95% CI. The heterogeneity 
of the effects across studies was assessed using the Chi square test 
and qualified by I 2 values, which represent the proportion of 
between-study variability that is attributable to heterogeneity 



rather than to chance; I 2 values of 25%, 50%, and 75% are 
referred to as low, moderate, and high heterogeneity effects, 
respectively. A fixed-effects model was used when there was no 
significant heterogeneity, whereas a random-effects model was 
used in other cases. Subgroup analyses were performed to 
investigate potential sources of heterogeneity by stratifying the 
different anti-TNFa agents, response criteria, duration of follow- 
up, and ethnic group. In addition, a sensitivity analysis was 
performed by the sequential omission of individual studies. 
Publication bias was estimated by funnel plots and an Egger's 
test, with p<0.05 considered as representing significant bias. All of 
the statistical analyses were performed using Review Manager 5.0 
software, and the results of the publication bias were estimated 
using ST ATA 12.0 software. 

Results 

Characteristics of the included studies 

Figure 1 describes the flowchart of our literature search. 
According to pre-determined criteria, fifteen articles, published 
between 2006 and 2012, were included [7,9-19,21-23] (A list of 
articles excluded and the reasons for exclusion were presented in 
Table S3 in File SI). As two described the same cohort [11,21], 
fourteen studies were ultimately included in this meta-analysis 
[7,9,10,12-19,21-23]. These studies comprised a total of 5561 RA 
patients: 5374 patients had their status of RF detected, and 1283 
patients had their status of anti-CCP antibody detected. All but 
one were prospective cohort studies. All of the participants in the 
14 studies had active RA, as evaluated by the number of swollen 
and tender joints, or DAS28 value, and the patients in all but 
one_ENREF_7 study had exhibited failure for at least one type of 
Disease-modifying Anti-Rheumatic Drug (DMARD) treatment 
[13], Almost all of the patients had established RA, with a mean 
disease duration from 5 to 14 years. The data extracted from the 
14 studies are summarized in Table 1. With regard to the anti- 
TNFa agents utilized, 10 studies involved infliximab, 8 involved 
adalimumab, 6 involved etanercept, and 1 involved golimumab. 

There was no study that involved patients receiving certolizu- 
mab. Only four studies investigated the status of RF, with only two 
investigating anti-CCP; the other eight studies measured both RF 
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Study or Subgroup 



Risk Ratio 
M-H. Random. 95% CI 



Risk Ratio 
M-H. Random. 95% CI 



Alexandras2009 

Bobbio2006 

BraunMoscovici2006 

Cuchacovich2008 

Lequerre2006 

Lilian2011 

Potter2009 

Vasilopoulos201 1 

Wijbrandts2008 

Wouter2008 



0.92 [0.60, 1.39] 
0.79 [0.63, 1.00] 
0.37 [0.20, 0.66] 
0.83 [0.54, 1.27] 
0.84 [0.59, 1.20] 
1.41 [0.82, 2.40] 
0.96 [0.87, 1.07] 
0.67 [0.53, 0.85] 
1.41 [0.94, 2.10] 
1.01 [0.85, 1.21] 

0.88 [0.76, 1.03] 



Total (95% CI) 
Total events 

Heterogeneity: Tau 2 = 0.03; Chi 2 = 27.32, df = 9 (P < 
Test for overall effect: Z = 1 .59 (P = 0.1 1 ) 



0.001); I 2 = 67% 



B 



Study or Subgroup 



Risk Ratio 
M-H. Random. 95% CI 



I 1 

0.2 0.5 1 2 

Favours experimental Favours control 

Risk Ratio 
M-H. Random. 95% CI 



1.7.1 EULAR criteria 

Alexandras2009 0.92 [0.60, 1.39] 

Bobbio2006 0.79 [0.63, 1.00] 

BraunMoscovici2006 0.37 [0.20, 0.66] 

Lequerre2006 0.84 [0.59, 1.20] 

Potter2009 0.96 [0.87, 1.07] 

Wouter2008 1.01 [0.85, 1.21] 

Subtotal (95% CI) 0.86 [0.73, 1 .02] 

Total events 

Heterogeneity: Tau 2 = 0.02; Chi 2 = 13.30, df = 5 (P = 0.02); I 2 = 62% 
Test for overall effect: Z = 1.76 (P = 0.08) 

1.7.2 DAS28 criteria 

Lilian2011 0.94 [0.77, 1.15] 

Vasilopoulos201 1 0.67 [0.53, 0.85] 

Wijbrandts2008 1.41 [0.94, 2.10] 

Subtotal (95% CI) 0.93 [0.65, 1 .34] 

Total events 

Heterogeneity: Tau 2 = 0.08; Chi 2 = 1 1 .32, df = 2 (P = 0.003); I 2 = 82% 
Test for overall effect: Z = 0.38 (P = 0.70) 



1.7.3 ACR20 

Cuchacovich2008 0.83 [0.54, 1.27] 

Lilian2011 1.41 [0.82, 2.40] 

Subtotal (95% CI) 1 .05 [0.62, 1 .79] 

Total events 

Heterogeneity: Tau 2 = 0.09; Chi 2 = 2.39, df = 1 (P = 0.12); I 2 
Test for overall effect: Z = 0.20 (P = 0.84) 



58% 



-t- 



-+- 



0.2 0.5 1 2 5 

Favours experimental Favours control 

Figure 3. Forest graphs of the meta-analysis of anti-CCP antibody status and response to anti-TNFa agents. The overall analysis of anti- 
CCP antibody status showed a pooled RR of 0.88 (95% CI: 0.76-1.03, p = 0.1 1 ) and an I 2 of 67%. Subgroup analyses on different response criteria 
revealed no significant differences. 
doi:10.1371/journal.pone.0089442.g003 
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Table 4. Subgroup meta-analysis of anti-CCP and RA patient response according to different anti-TNFa agents, follow-up periods, 
response criteria, and ethnic groups. 













Test of association 




Test of heterogeneity 




Subgroup 


Population 


No. of 
studies 


No. of 
patients 


RR 


95%CI 


P value 


Model 


P value 


I 2 




overall 


10 


1283 


0.88 


0.76-1.03 


0.11 


R 


<0.01 


67% 


Anti-TNFa agent* 


infliximab 


5 


345 


0.79 


0.55-1.13 


0.19 


R 


<0.01 


78% 




adalimumab 


3 


297 


1.01 


0.88-1.17 


0.86 


F 


0.69 


0% 


Follow up period 


5:6 months 


7 


1074 


0.90 


0.80-1 .01 


0.06 


R 


0.09 


46% 




<6 months 


3 


209 


0.78 


0.40-1 .53 


0.47 


R 


<0.01 


86% 


Ethnic group 


European 


7 


1144 


0.90 


0.79-1.04 


0.17 


R 


0.02 


60% 




South American 


2 


109 


1.01 


0.79-1.30 


0.91 


F 


0.38 


0 




Asian 


1 


30 


0.37 


0.20-0.66 


<0.05 


NA 


NA 


NA 



F, fixed-effects model; R, random-effects model; NA, not applicable. 

*The number of studies and number of patients receiving individual anti-TNFa agent in this table are those that we could identified after literature review and 

contacting the authors. 
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and anti-CCP. The follow-up duration of these 14 studies ranged 
from 14 weeks to 48 weeks. 

Study quality 

Our assessment was performed on the basis of Hayden's criteria 
for prognostic studies [20]. The results revealed that only one of 
the included studies was considered to be of low quality, with a 



score of 7; the other thirteen studies were considered as having 
moderate or high quality (1 1 moderate and 2 high) (Table 2). 

Association evaluation of RF/anti-CCP antibody status 
and patient response to anti-TNFa agents 

Twelve of the included studies, with a total of 5374 patients, 
examined the predictive effect of the RF status for patient response 
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Figure 4. Overall analysis of publication bias on the effect of RF status on the response to anti-TNFa treatment. Egger's linear 
regression test was performed to quantify publication bias. The p values of the RF status analysis were 0.777. The funnel plot showed no significant 
evidence of asymmetry. 
doi:1 0.1 371 /journal.pone.0089442.g004 
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Egger's test 
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Figure 5. Overall analysis of publication bias on the effect of anti-CCP antibody status on the response to anti-TNFa treatment. 

Egger's linear regression test was performed to quantify publication bias. The p values of the anti-CCP antibody status analysis were 0.422. The funnel 
plot showed no significant evidence of asymmetry. 
doi:1 0.1 371 /journal.pone.0089442.g005 



to anti-TNFa agents. The overall analysis revealed that the pooled 
RR was 0.98 (95% CI: 0.91-1.05, p = 0.54) and that I 2 was 43% 
(Figure 2a), suggesting that the RF status was not associated with a 
patient's response to anti-TNFa treatment, with a moderate 
heterogeneity observed. Subgroup analyses on the different anti- 
TNFa agents, follow-up periods, response criteria, and ethnic 
groups were performed to investigate sources of this heterogeneity, 
though the results revealed no significant differences (Figure 2b, 
Table 3). 

Ten studies, with a total of 1283 RA patients, were analyzed for 
an association of anti-CCP antibody status and patient response to 
anti-TNFa treatment. The overall meta-analysis showed no 
association between the status of anti-CCP antibody and a 
patient's response to anti-TNFa treatment, with a pooled RR of 
0.88 (95% CI: 0.76-1.03, p = 0.11) and an I 2 of 67% (Figure 3a). 
When we performed a subgroup analysis to explore potential 
heterogeneity, we did not identify any association by stratifying the 
studies using either different anti-TNFa agents, follow-up periods, 
response criteria, or ethnic groups (Figure 3b, Table 4). 

Sensitivity and publication bias 

Sensitivity analyses were performed based on the results of the 
bias assessments, Hayden's criteria of quality assessment for 
prognostic studies, by the sequential omission of individual studies, 
and the result revealed that the significance estimate of the overall 
pooled RR was not influenced by omitting any single study. The 
publication bias was evaluated by a funnel plot, which showed no 
significant evidence of asymmetry (Figures 4-5). We also 
performed an Egger's linear regression test to quantify the 
publication bias. The p values of the RF and anti-CCP antibody 



status analyses were 0.777 and 0.422 (Figures 4-5), and all other 
Egger's test p values in this meta-analysis were >0.05, suggesting 
no significant bias in the analysis. 

Discussion 

Previous single studies have suggested that the presence of RF 
and anti-CCP antibody is associated with a reduced response to 
anti-TNFa treatment [7-14];_ENREF_9 however, a number of 
studies have reported a contrasting conclusion [15-19]. Thus, to 
investigate whether the status of RF or anti-CCP antibody in RA 
patients has predictive value for a clinical response to anti-TNFa 
treatment, we performed a meta-analysis by a comprehensive 
literature search and inclusion of all of the available qualifying 
studies. The results of our meta-analysis demonstrated that there 
was no association between the status of RF and patient response 
to anti-TNFa treatment. Furthermore, the status of anti-CCP 
antibody was also not associated with patient response, with a 
pooled RR of 0.88 and 95% CI of 0.76-1.03 (p = 0.11), indicating 
a tendency of association between the absence of anti-CCP 
antibodies and a patient's response to anti-TNFa therapy. 
Therefore the cumulative results did not support or completely 
reject the association between the presence of RF and/ or anti- 
CCP and the effect of treatment with TNFa inhibitors. Anti-CCP 
is an independent predictor of radiological damage and progres- 
sion in RA patients [24,25], and anti-CCP-positive patients have a 
relatively worse prognosis. As a result, patients with a different 
anti-CCP antibody status may have a different baseline of the 
disease, which could influence their response to anti-TNFa agents. 
To explore this influence, analyses should be performed to 
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investigate potential differences in the baseline characteristics 
between groups stratified by the status of anti-CCP antibody. 
Alternatively, only analyses based on studies with placebo groups 
are powered to definitively identify such an association. Unfortu- 
nately, such information was not available for the included studies. 

This meta-analysis does have some limitations. First, only 14 
studies were included in this meta-analysis, with 12 studies 
reporting the status of RF and 1 0 studies reporting the status of 
anti-CCP antibody. Nonetheless, after a comprehensive literature 
search covering five databases and the selection of eligible studies 
by three different investigator strictiy according to the inclusion 
criteria, most of the included studies were assessed with moderate- 
to-high quality. Therefore, it is reasonable to obtain a conclusion 
based on the studies included in this meta-analysis. However, as 
potential heterogeneity exists between studies after subgroup 
analyses were performed, the result of our analyses should be 
interpreted with caution, and our conclusions should be updated 
using new, larger studies in the future. Secondly, all of the patients 
in the studies included in our meta-analysis had long-standing 
disease, with a mean value of disease duration ranging from 5 to 
14 years, and had previously shown failure with at least one type of 
DMARD. Indeed, patients with more severe chronic disease as a 
result of irreversible joint damage may be less likely to respond to 
anti-TNFot treatment [17]. Thirdly, most of the included studies 
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